Melting of Flux Lines in an Alternating Parallel Current 
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We use a Langevin equation to examine the dynamics and fluctuations of a flux line (FL) in 
the presence of an alternating longitudinal current J\\(uj). The magnus and dissipative forces are 
equated to those resulting from line tension, confinement in a harmonic cage by neighboring FLs, 
parallel current, and noise. The resulting mean-square FL fluctuations are calculated exactly, and 
a Lindemann criterion is then used to obtain a nonequilibrium 'phase diagram' as a function of the 
magnitude and frequency of Ju (ui). For zero frequency, the melting temperature of the mixed phase 
(a lattice, or the putative "Bose" or "Bragg Glass") vanishes at a limiting current. However, for 
any finite frequency, there is a non-zero melting temperature. 
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I. INTRODUCTION 



The phases of flux lines (FLs) in high temperature su- 
perconductors are the subject of many current experi- 
mental and theoretical investigations [Qj2|. In the classi- 
cal (Abrikosov) picture of type-II superconductors, FLs 
penetrate a clean material for fields H > H c i , to form a 
triangular solid lattice. This mean field phase diagram is 
modified by inclusion of thermal fluctuations, or in the 
presence of various forms of disorder (such as point or 
columnar defects), which play an important role in high- 
T c materials (|||. The essentially classical problem of 
understanding the equilibrium and non-equilibrium prop- 
erties of the resulting vortex matter composed of mean- 
dering FLs, poses a myriad of interesting challenges [Q-Q . 

There is now growing evidence that the mixed phase of 
high-T c materials is a "Bose Glass" (BoG) [^P], in which 
the FLs are pinned by such extended (or correlated) im- 
perfections as grain boundaries or columnar pins. Since 
the FLs are pinned to these large defects, the BoG has an 
infinite tilt modulus C44. There is also a cusp- like singu- 
larity in its melting temperature in a field H± , perpendic- 
ular to the extended defects || . On the other hand, if the 
defects are at uncorrelated points (such as oxygen impu- 
rities) [Q, there may exist a "Bragg Glass" (BrG) phase 
at low magnetic fields and temperatures || . The disloca- 
tion free "Bragg Glass" has a kind of topological order. 
Experiments on BSCCO || have measured the effects of 
both point and columnar imperfections, showing that a 
melting transition exists in the presence of either type of 
disorder. For technical applications, it may be useful to 
control the melting temperature of the mixed phase (BoG 
or BrG) by application of external parameters. One pos- 
sibility is an alternating longitudinal current J»(u>), as in 
some recent experiments JToj ] . 

The response of FLs to Jy was studied extensively even 
prior to the high-T c era H. Since the Lorentz force 



F = (f)o J x B is absent for J || B, the straight FL is 
stable. Any fluctuation of the FL, however, results in a 
force that causes helical instabilities at long wavelengths 
B . A lattice of FLs is similarly unstable to an infinites- 
imal Ju , if it penetrates the whole sample. As shown by 
Brandt |11|, for a finite penetration depth of current, a 
finite Ju ~ e 1 / 4 , is needed to cause instability, where e 
is related to the shear modulus cqq. The many studies 
of this subject point to the importance of surface effects 
[^],fLlJ, and the bewildering behavior of the unstable FLs. 

Brandt has suggested that FLs can also be stabi- 
lized for infinitesimal Jm due to pinning by disorder. How- 
ever, this study predates the introduction of the BoG and 
BrG phases, which are more ordered than a flux liquid. 
We may thus enquire if these phases are stable with re- 
spect to J|| . More generally, what is the behavior of the 
flux matter in response to an alternating current J||(w)? 
In what follows, we shall assume a current Ju (u>) that 
uniformly penetrates the whole sample, and calculate a 
'melting temperature' (T m ) based on this assumption. 
Self-consistency of the treatment thus requires penetra- 
tion of Ju through the whole sample. Whether or not 
the longitudinal penetration depth Ay (u>) is finite at T m 
appears to be unsettled. This will be discussed in more 
detail in the next sections. If it does happen that A|| is fi- 
nite at the transition, our result becomes applicable only 
to thin samples with thickness of a few An, or alterna- 
tively, only close to the surface of the sample. Nonethe- 
less, such a melting at the surface will undoubtedly effect 
the transport properties of the superconductor, and may 
well trigger melting in the bulk of the material. 

We use the cage model to examine the fluctuations of 
one representative FL in the confining potential provided 
by its neighbors. The collective interactions of many FLs 
is thus approximated by a simple harmonic potential. 
The cage model provides a good description for the BoG 
phase in the limit of B > 0o/A^ as argued in Ref. ||, 
and is also applicable to the BrG, as discussed by Er- 
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tas and Nelson @. Basically, there is evidence of some 
reasonably long, but possibly finite range, translational 
order in the flux phase of high temperature materials, 
despite the presence of disorder (^J^j. The cage model 
is a valuable tool for estimating the limits of such order, 
and how it is destroyed by fluctuations. More recently, 
the cage model was used by Goldschmidt |l2|, in con- 
junction with a replica variational method, to compare 
and contrast the effect of point and columnar disorders 
on the melting temperature. 

So far the cage model has only been used |3|,|| to study 
melting due to thermal fluctuations or impurities under 
equilibrium conditions. The Lindemann criterion, how- 
ever, has been applied to study the non- equilibrium melt- 
ing dynamics of the vortex lattice moving in a random 
medium due to a transverse current Jl3[ . In the same 
spirit, we use the Lindemann criterion, in conjunction 
with the cage model, to examine the destruction of order 
induced by an alternating parallel current. Inclusion of a 
time dependent current necessarily implies the need to go 
beyond the standard treatments of a Hamiltonian . 
Instead, we work directly with a Langevin equation for 
the motion of a FL, enabling us to estimate fluctuations 
in the presence of J\\(lo). 

The rest of this paper is organized as follows. In Sec. II, 
we introduce the cage model and construct an equation 
of motion for the FL in the presence of J\\(u)- Using a 
Lindemann criterion, we then estimate a melting tem- 
perature T m , and examine its dependence on the ampli- 
tude and frequency of J|| . Application of this result to 
BoG and BrG are discussed in Sees. Ill and IV respec- 
tively. The conclusions appear in sec. V, and details of 
the steady state distribution for lu = are placed in the 
Appendix. 



II. THE CAGE MODEL 

In the spirit of the Einstein model for vibrations of a 
solid, the "cage model" approximates the system of many 
interacting FLs by that of a single FL confined by its 
neighbors to a harmonic potential kr 2 /2. The spring con- 
stant fc, is approximated in the dense limit of B > H Cl , 
by k w £o/a 2 ; where £o = (0o/47rA) 2 , and (f>o = hc/2e. 
The spacing of flux lines is indicated by a, while A is the 
penetration depth |IJ . The energy of the FL is then given 
by the Hamiltonian 



H= dz 
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(1) 



where r (z) is a two-component vector describing the FL 
configuration (i || B). In the addition to the harmonic 
potential, there is an effective line tension e\ = £0/7 , 
where 7 2 = m z /m± ^> 1 is the mass anisotropy. 

In the Langevin approach, the dynamics of the field 
r(z,t) is governed by the equation of motion 
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r = -[i-^ + ff(z,t) = p (eid z 2 r- kf) +ff(z,t), (2) 

where fi is the mobility of the FL, estimated by the 
Bardeen-Stephen |l5|] expression fi ~ p n /(cf)oH C2 ), where 
p n the normal state resistivity of the material. The ran- 
dom noise if has zero mean, and is uncorrelated at differ- 
ent points and times, such that 

(ri a (z, t)rjp(z', t')) = 2D5 af3 5(z - z')S(t - t'). (3) 

To ensure that in the steady state the fluctuations are 
governed by the Boltzmann weight of the Hamiltonian Ti. 
at a temperature T, we must satisfy the relation 



D = fiksT. 



(4) 



The actual dynamics of a FL can be more complicated. 
Even if restricted to local terms, the symmetries of the 
system allow the inclusion of a Magnus force, generalizing 
Eq. (0) to 



r + az x r = — ^-r^r + i](z, t), 
or 



(5) 



where fi is again the mobility. The parameter a is related 
to the Hall angle in the normal phase by a = tan 9h [0 , 
and expresses the tendency of the FL to drift perpendic- 
ular to the applied force. In the Appendix, we prove that 
even for this more general equation, the requirement of 
thermal equilibrium is still satisfied if Eq. (|) holds. 

A time dependent current J(i), destroys thermal equi- 
librium and can not be placed in the Hamiltonian. Such 
a current generates a Lorentz force F = </>o J X t, where t 
is the unit tangent vector. In the case of a current J || z, 
the (linearized) force is 4>oJ z (t)z x d z f. We shall assume 
that close to equilibrium this force can be simply added 
to the time evolution in Eq. (||) . Thus we shall study the 
equation of motion 

r + az x f— fi [sid z 2 r — kf+ 4> J z (t)z x d z f ] + rf, (6) 

while keeping the fluctuations of the noise related to am- 
bient temperature by Eqs. (||) and (Q). 

By decomposing the vector f(z,t) into its two compo- 
nents r\ and T2, and using the transformation R(z,t) — 
ri(z,<) + ir2(z,t\ Eq. can be written as 

(1 + ia) R = n [e!d z 2 R -kR + i(j) J z (t)d z R] + fj, (7) 

where 77 = 771 + ir/2. After transforming to Fourier space, 
we obtain 

(1 + ia) R + n[e x q 2 + k + 4> J z (t)q] R = rj(q,t), (8) 
whose solution is given by 
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(9) 



exp ^~ l ^ ia dT [ £ i1 2 + k + ^Jz{r)q] 
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where we have assumed that R(q,0) = 0. 
From Eq. (§), 

we obtain that the Fourier components 

of noise satisfy 

(fj(q,t)f}*(q',t , ))=4D(2w)S(q-q , )5(t-t'). (10) 
Using Eqs. @ and (10) we find, 



i\R( q ,tr) = fdf^x 

Jo 1 + a 1 



(11) 



exp i 1 + 2 ^ 2 jf dr [eiq 2 + k + <f>oJ z {T)q] 



where we have set 27r<5(0) = L z , for a sample of length 
L z in the z-direction. 

We shall henceforth consider an alternating current, 
Jz — JzCOs(wi), for which the integral in the exponent 
can be easily performed to yield 



(\R(d,t)\ 2 ) 



dt' 
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(12) 



/ 9 , \ / is , sin(cji) — sin(wt') 
(e iq 2 + k)(t- t') + <j> j z q^^ 5: — L 



Note that the mean-square fluctuations have an oscilla- 
tory dependence on the time i, even in the long-time 
limit. However, to calculate the fluctuation-induced 
melting temperature we have to consider the behavior 
over the entire cycle. In particular, we shall focus on the 
maximal fluctuations which occur when tut is a multiple 
of 7r |16[ ]. After setting t' = t — t', we obtain in the 
long-time limit 



(|iW<z)l 2 ) = T 



dr x 



(13) 



exp 



-2/zt' 

TT^ 2 



, ^OjzQ . i i 

E\q + k ± — sin(ur ) 



by 



The extent of FL fluctuations in position-space is given 



I I f{r l {zY + r 2 {zf) = ± I ^(\R(q)\ 2 
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(14) 



Following simple changes of variables, the maximal ex- 
tent of fluctuations is obtained from Eq. (|l3| ) as 



I = Ir 



where I c 



f°° dx 




L vi exp i 






r 



! , Jz\ 2 ( sin(wa;/w c ) x ' 



lox/lo c 



(15) 



D/fiy/wkei w D4> H C2 a^f / '{y/np n eo), j c 



^kei/(j) Q 2 w 2e /(a7? !, o), and cj c = 2/ifc/(l + a 2 ) = 
2p n Eo/{a 2 <p H C2 {l + a 2 )). Using Eq. (§), we find that 7 C 



depends on the temperature as I c = ksT/y/nkei, but is 
independent of a, while w c does depend on a. 

According to the Lindemann criterion, melting occurs 
when sj (f 2 ) ps clcl, where cl is a constant of the order 
of unity. Using this criterion, we can find the dependence 
of T m on j z and lo. For ui = the integral in Eq. (|l5| ) 
is given by J = /cV^/fl = Uz/jc) 2 } for j 2 < j c , and di- 
verges for j z > j c . Therefore, the melting temperature 
is 



m. in 
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jc 



1/2 



for 



Jz < Jc 
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as 



where ksT^ = 

parameters for BSCCO, 7 = 50, 

100(K/A), a = 225A(B w 4T), p n = 8~x 10" 3 (f7 - cm), 
and a = 1.75, leads to j c s» 40(^4/m 2 ), and w c 17kHz. 
Finally, the choice of c L = 0.2, leads to T,° » 18°K, 
which is the estimate for the melting temperature in the 
absence of a current J z . (Similarly, the parameters 7 ps 7 
and e « 50 appropriate to YBCO, yield T° t w 64°^.) 

The divergence of the integral (for j z > j c ) is removed 
at any finite frequency u>, leading to a non-zero melting 
temperature. The behavior of T m for different values of j z 
and ui (obtained from numerical integration) is shown in 
Fig. 1. Due to the onset of helical instabilities for ui = 0, 
the melting temperature vanishes for j > j c . However, 
at any finite frequency the system retains its stability 
for all currents, although fluctuations are enhanced. The 
cage model provides an estimate for the destruction of 
translational order by these fluctuations. 
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Fig. 1 Dependence of the melting temperature on j z and 
u>. The solid line is for u> = 0, (a) for ui/lu c = 0.5, and (b) for 
ui/uJc = 1. 



III. THE BOSE GLASS 

To what extent is the previous discussion of a single FL 
in a cage applicable to the disordering of the flux phases 
of high temperature superconductors? The susceptibility 
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of FLs to pinning by the many sources of impurities in 
these materials destroys the long-range translational or- 
der of a FL lattice || , replacing it by one of a number of 
potential glassy phases. These glass phases have vanish- 
ing linear resistivity because of the presence of diverging 
energy barriers at low current density || . From the the- 
oretical perspective, the best understood type of glass 
phase results from correlated disorder, such as parallel 
columnar pins or grain boundaries. Because of a map- 
ping H from this anisotropic system to superfluid bosons 
on a disordered substrate, the resulting phase is called a 
Bose Glass (BoG). Due to its anisotropy, the BoG phase 
has a distinct signature in its conductivity as a function 
of orientation. While the majority of transport experi- 
ments are carried out with currents perpendicular to the 
magnetic field, recent experiments [ [To| measuring elec- 
trical transport parallel to the field direction, appear to 
lend support to the BoG theory. 

As shown by Nelson and Vinokur H , the BoG phase is 
distinct from a liquid of FLs, and characterized by a rea- 
sonably long translational order despite the presence of 
columnar defects. The destruction of this order upon in- 
creasing temperature (or magnetic field) is reminiscent 
of melting, occuring at an irreversibility line which is 
typically at higher temperatures than that of pure lat- 
tice melting. It has been argued in Refs. pUl7| that the 
cage model provides a good description for this transition 
when the magnetic field is greater than — <po/X^_, 
which is about IT for BSCCO. The results of Fig. 1 
should thus apply to the BoG phase in this region, es- 
tablishing the limits of its stability in the presence of a 
longitudinal current. 

For j z — 0, Eq. ( |l6| ) reproduces the equilibrium melt- 
ing temperature of T rn = T® n = ac^eo/^f, as in previous 
studies H- However, applying the cage model at finite 
J|| (lu) assumes that the current uniformly penetrates the 
whole sample. Self-consistency then requires that the 
current penetration depth should diverge at the melting 
temperature; an issue which is not quite clear at present. 
Ref. H claims that when T — ► T^ oG , the penetration 
length goes to infinity. However, Ref. |l8| finds that (for 
u> — > 0) the penetration length diverges as T — > Tg oG for 
a perpendicular current, but remains finite for parallel 
currents. (Note that unlike usual melting, the disorder- 
ing of the BoG phase occurs through a continuous phase 
transition.) The AC current also penetrates the sample 
at the scale of a few A for the relatively (to the gap) low 
frequencies close to ui c . If the current does not fully pen- 
etrate the sample, the bulk of the material may remain in 
the BoG phase, while the surface layer melts. The melt- 
ing of the FLs at the surface will however profoundly 
effect the transport properties of the material, even if no 
structural change is detected in the bulk. 



IV. THE BRAGG GLASS 

The nature of glassy phases in the case of isotropic dis- 
order (e.g. point defects due to oxygen vacancies) is more 
complex. One proposal is for a "gauge glass" which as- 
sumes complete destruction of the lattice order. Another 
possibility is a dislocation-free Bragg Glass (BrG) at low 
temperatures and magnetic fields. The BrG maintains 
the topological order of a lattice, and the corresponding 
Bragg peaks in scattering experiments M . Ref. || argues 
that due to its similarity to a solid, interactions of a FL 
with its nearest neighbors in the BrG are well captured 
by the cage model. Therefore, one can again apply our 
results to the BrG phase, and the phase diagram of Fig. 1 
can also be used to establish the limits of stability of this 
phase in the presence of a current J\\(lu). 

In the limit of j z = 0, the equilibrium melting tem- 
perature of T m = 2a£ocf j /7 from Eq. jl6[ ) is consistent 
with the result of T m = 4a 3 ^0^044 c 2 L in Ref. . The 
latter is obtained using replica and variational methods, 
and coincides with the former if we set cqq = Eq /4a 2 and 
C44 = Eq/ J 2 a 2 (single vortex contributions). The pro- 
visos presented before regarding the penetration of the 
alternating current into the bulk also apply to this case. 
However, the melting of the BrG is thought to be first 
order, excluding the possibility of diverging penetration 
depths at the transition. 



V. CONCLUSIONS 

In this paper we examine the response of a single FL 
in a harmonic cage to a (destabilizing) parallel alternat- 
ing current J\\(lu). A Langevin equation is used to track 
fluctuations of the FL in the presence of both Jn (lu) and 
thermal noise at a temperature T. The harmonic cage is 
frequently used to mimic the confinement of a FL by its 
neighbors in ordered flux phases. The equilibrium melt- 
ing temperature (for the destruction of such order) is then 
estimated by the Lindemann criterion which equates the 
extent of fluctuations to a fraction of the FL spacing. 
While the system is no longer at equilibrium in the pres- 
ence of an alternating current, we can still ask if the FL 
order is maintained in the steady state. In regions where 
the current penetrates the sample uniformly, we again 
apply the Lindemann criterion to answer this question. 

Our primary results are summarized by Fig. 1: (1) 
At zero temperature and frequency, the FLs are stable 
up to a current j c which is related to the strength of 
the confining potential. (2) A uniform current (lu = 0) 
decreases the energy cost of long wave-length deforma- 
tions, resulting in larger FL fluctuations in response to 
thermal noise. The melting temperature T m (J||) is thus 
depressed, eventually going to zero for Jm = j c . (3) The 
diverging response of the FL is removed at any finite 
frequency lu. Thus while the transition temperature is 
reduced by the alternating current, it remains finite. A 
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substantial reduction of the transition temperature oc- 
curs only at frequencies less than a characteristic value 
of uj c related to the FL mobility (and hence to the re- 
sistivity in the normal state). For BSCCO we estimate 
j c sa 4Q(A/m 2 ), and uj c ss 17kHz. 

While the Lindemann criterion is usually applied to a 
solid phase, it should be equally applicable to any of the 
glassy phases which have a substantial degree of trans- 
lational order. Our results are thus also relevant to the 
dislocation free BrG, as well as to the BoG in the dense 
limit. While the pinning of FLs by disorder stabilizes the 
glassy phases to the current, we note that the ordered 
lattice in the absence of impurities is in fact unstable for 
even infinitesimal J\\{lu = 0). It would be interesting to 
examine the stability of a collection of FLs in the pres- 
ence of both impurities and parallel current by extending 
previous replica variational calculations Even the 

fluctuations of a single FL in the absence of a cage to the 
combined effects of current and randomness is an inter- 
esting topic which can be potentially studies by computer 
simulations. Another potential extension is to study the 
depinning transition of the FL at zero temperature by 
methods such as in Refs. fl[l9|. 

Finally we note that a related situation arises in the 
context of chiral polymer crystals . Long polymers in 
dense solution often crystallize into a hexagonal colum- 
nar phase. When the polymers are chiral, a bias towards 
cholesteric twist competes with braiding along an average 
direction. If the chirality is weak, a defect free hexago- 
nal columnar phase persists, as in the Meissner phase. 
When the chirality is strong, screw dislocations prolifer- 
ate, leading to either a tilt grain boundary phase or a 
new state with twisted bond order. This new phase may 
be similar to a flux lattice in the presence of a parallel 
current |21|. 
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M. Kohandel acknowledges support from the Institute 
for Advanced Studies in Basic Sciences, Gava Zang, Zan- 
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Appendix 

Despite the differences in Eq. (^) from a standard 
Langevin equation, the calculated averages of (r 2 ) for 
u> = are those from the Boltzmann weight of a Hamil- 
tonian 



(6i 



SH 

OTi 



(18) 



n 



dz 



»o3 



H 7T~ ( r 2<9 z ri - nd z r 2 ) 



(17) 



First, we note that a uniform current can be absorbed 
into an effective Hamiltonian as also found in Ref. ||. 
Second, it can be shown quite generally that a Langevin 
equation of the form 



where is the anti-symmetric tensor, has a station- 
ary state in the form of a Boltzmann weight po oc 
exp [— H/ltsT], provided that the correlations of the 
noise satisfy the condition in Eq. (||). For other ex- 
amples of fluctuation-dissipation conditions in modified 
Langevin equations, see Ref. |22]. 
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